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ABSTRACT 

Feedback in massive galaxies generally involves quenching of star formation, a favored candidate 
being outflows from a central supermassive black hole. At high redshifts however, explanation of the 
huge rates of star formation often found in galaxies containing AGN may require a more vigorous 
mode of star formation than attainable by simply enriching the gas content of galaxies in the usual 
gravitationally-driven mode that is associated with the nearby Universe. Using hydro dynamical sim¬ 
ulations, we demonstrate that AGN-pressure-driven star formation potentially provides the positive 
feedback that may be required to generate the accelerated star formation rates observed in the distant 
Universe. 

Subject headings: galaxies: formation — galaxies: active — methods: numerical 


1. INTRODUCTION 

The remarkable universality of the Schmidt-Kennicutt 
star formation law, ranging from global fi ts to star¬ 
formi ng galaxies in the nearby Universe (Kennicutt 
|1998|), local fits to star-for ming compiexes within galax¬ 
ies (iKenni cntt et ai.||2QQ7|), an d to star-forming galaxies 
to z ^ 2, (|Genzel et al.||201Up inspires considerable con¬ 
fidence in the theory of gravitationai instabiiity-driven 
star formation in gaiactic disks. There are exceptions, 
most notabiy in molecui ar complexes with ano maiousiy 
low star formation rates (Rathborne et al. 2014), but sim¬ 
ple and plausible additions to the usual density threshold 
criterion for star formation, most notably by incorporat¬ 
ing turbulence, may go far towards resolving these i ssues, 
as demonstrated both theoretica lly (for a review see Krit- 
suk, Ustyugov & Norman||2011 ) and pheno menologically 
in we ll resolved examples such as NGC 253 (Leroy et al. 


In the high redshift Universe, the accumulation of re¬ 
cent data on remarkably hig h star formation rates poses 
a fascinating challenge, (e.g. Drouart et al.||2Q14 Picon- 
celli et al. 2015 Rodighiero et al. 2015). Is it simply a 
question of turning up the gas traction or is a new mech¬ 
anism at work for inducing more efficient star formation? 

In this Letter, we reinforce the case for the latter, more 
radical view by simulating the evolution of a fully self- 
consistent, gas-rich star-forming disk galaxy that is sub¬ 
ject to the overpressuring influence of a vigorous outburst 
from its central Active Galactic Nuclei (AGN). Star for¬ 
mation is a complex interplay between gas supply, multi¬ 
phase interstellar medium (ISM), cloud collapse, and gas 
ejection from the disk. We will consider two cases, cor¬ 
responding to gas-poor (gas fraction 10%) and gas-rich 
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(gas fraction 50%) systems at z ^ 2. 

The case for inducing star for mation via AGN activ¬ 
ity has been ma de analytically (Sil k fc Rees 199^ |Silk 
fc Norman|2009 ) and in s imulations ([Gaibler et al.|2012[ 


Ishibashi ^ habiai^|2012ljWagner, Bicknell 


memura 


[Zubovas et al.||2013j) with varying degrees of astro 

physical reality, birnulations demonstrate that whether 
the AGN activity is jet or wind-induced is irrelevant: 
after a few kpc, both inputs are indistinguishable. Hith¬ 
erto, however, star formation in the multiphase ISM has 
not been followed in adequate detail because of the lack 
of numerical resolution. 

For the first time, we study the effects of pressurization 
of the disc by performing simulations in which we fully 
include self-gravity of the multiphase ISM and thereby 
trace the evolution of the star formation rate (SFR) as 
well as that of the gas content of the system. 

2. SIMULATION SET-UP 

To study the efffect of an external pressure on a galaxy, 
we have performed a four isolated disc galaxy simulations 
with two different initial gas fractions of 10% (hereafter, 
gasLow) and 50% (hereafter, gasHigh). We allow the 
galaxies, of one-tenth the total mass of the Milky Way, 
to initially adiabatically relax to an equilibrium configu¬ 
ration (with a reasonable disc thickness) over the rotation 
time of the disc at its half-mass radius. After this first 
relaxation phase, we turn on the external pressure, gas 
cooling, star formation, and also feedback from super¬ 
novae (SNe), as will be described below. 

In the initial conditions, the dark matter (DM) parti¬ 
cles a re sampled with an NFW (Navarro, Frenk & White 
1997|) density profile and a concentra tion parameter oT 
c = 1 0 using the method introduced by |Springel fc Her~ 


quist (2005). For the DM particles, the virial veloc 
ity is set to be i; 2 oo = 70kms“^, which corresponds to 
a virial radius of i ?200 ~ 96 kpc and a virial mass of 
M 200 ~ 1.1 X 10^^ Mq. We use 10^ DM particles with a 
mass resolution of 1.21 x 10^ Mq to sample the DM halo. 
These stellar disk and bulge were initially sampled with 
5.625 X 10^ particles of which 6.25 x 10^ were used to 
sample the bulge. The stellar particles are distributed in 
an exponential disc with a scale length of 3.44 kpc and 
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scale height 0.2 kp c, and a spherical, non rotating bulge 
with a Hernquist (1990) profile of scale radius 0.2 kpc. 

The simulations are run with the R AMSES adaptive 
mesh refinement code (Teyssier 2002). The box size 
is 655 kpc with a coarse level of 7, and a maximum 
level of 14 corresponding to a maximum resolution of 
Ax = 40 pc. The refinement is triggered with a quasi- 
Lagrangian criterion: if the gas mass within a cell is 
larger than 8x10^ Mq or if more than 8 dark matter 
particles are within the cell a new refinement level is 
triggered. 

After relaxation, the origin of time was reset to zero 
and the base simulations were run further in time with 
an enhanced and uniform pressure outside the disc {pres¬ 
sure simulations) for another 0.42 Gyr. This pressure 
enhancement is applied instantaneously at t = 0 for a 
value of 3Pmax (hereafter pa3) outside the sphere of ra¬ 
dius ri = 12 kpc. Pmax is the maximum radially averaged 
pressure in the disk at t = 0 (reached in the central few 
cells), with Pmax — 9.8 x 10“^^ Pa for the gasLow simula¬ 
tion set and 4.7 x 10“^^ Pa for the gasHigh simulation 
set. 

Sub -grid models for cooling Sutherland fc Dopita| 
(1993) and star formation, as well as SN feedback, were 
used in the simulations. Gas is turned into star parti¬ 
cles in dense cold regions of gas density rigas > no = 
14Hcm“^ by drawing a probability from the Schmidt 
law p* = O.Olpgas/te to form a star with a stellar mas s 
of = tiqAx^ 2 X 10^ Mq (Rasera & Teyssier||2006 ). 
The gas pressure and density are evolved using the Euler 
equations, with an equation of state for a mono-atomic 
gas with 7 = 5/3. In order to prevent catastrophic and 
artificial collapse of the self-gravitating gas, we use a 
polytropic equation of state T = To(ngas/no)^“^ to arti¬ 
ficially enhance the gas temperature in high gas density 
regions (ngas > no). Here k. = 2 the polytropic index, 
and To = 270 K, chosen to resolve the Jean s length with 
minimum 4 cells (jDubois & Teyssier||2008|). We account 
for the mass and energy release from type II SNe. The 
energy injection, which is purely thermal, corresponds to 
^SN = ^SN (ui^/Mq) 10^° erg, where psN = 0.2 is the 
mass fraction of stars going into SNe. We also return an 
amount rjsN'm^ back into the gas for each SN explosion 
which occurs 10 Myr after the birth of the star particle. 
To avoid excessive cooling of the gas due to our inability 
to capture the different phases of the SN bubble expan- 
sio n, we used the delay ed cooling approach introduced 
by Teyssier et al. (2013). 


3. RESULTS 

3.1. Dise fragmentation 

The application of external pressure at t = 0 leads to 
fragmentation of the gaseous galaxy discs, i.e. acceler¬ 
ated clump formation. This can be seen in Fig. for the 
gasLow and in Fig. for the gasHigh simulations. How¬ 
ever, the gasHigh case shows more gas between clumps in 
the enhanced pressure run, as well as more gas ejection 
from the disk. 

Since the star formation recipe depends on the local 
gas density, we expect an enhanced star formation when 
more clumps are formed (gas gets more concentrated), 
assuming that the clumps have sufficient mass. There¬ 
fore, if external pressure leads to increased fragmenta¬ 


tion and hence increased clump formation, we expect 
the star formation to be positively enhanced when exter¬ 
nal pressure is applied on the galaxy. We first consider 
the fragmentation by counting the high-density clumps. 
To find the clumps, we run the clump finder described 


Bleuler fc Teyssi^ (2014), which identifies all peaks 


and their highest saddle points above a given threshold 
(%as > 20Hcm“^). A clump is recognized as an individ¬ 


ual one when the peak-to-saddle ratio is above 1.5, oth¬ 
erwise the peak is merged with the neighbor with which 
it shares the highest saddle point. 

The top panel of Fig. shows the number of clumps 
as a function of time. One can see that generally more 
clumps are formed when external pressure is applied. In 
the gasLow simulations, the number of clumps starts ris¬ 
ing soon after this pressure is applied. In the absence 
of external pressure, the number of clumps is unity (the 
disc as a whole corresponds to a clump). This value is 
only reached at late times (after 300 Myr) once the gas 
has sufficiently collapsed to reach the clump gas density 
threshold. 

On the other hand, in the gasHigh simulations, clump 
formation sets in at earlier times, even without any forc¬ 
ing by external pressure. This diff e rent b ehavior is re¬ 
lated to differences in the Toomre (1964) Q parameter 
in the galaxies of the two simulations, such that the 
gaseous disc is stable in the gasLow case {{Q) = 3.29 > 1 
measured at t = 0), but unstable in the gasHigh case 
{{Q) = 0.72 < 1). This demonstrates that the fragmen¬ 
tation of the disc can be driven by the forcing of an ex¬ 
ternal pressure, even though the disc is initially Toomre 
stable. For the gasHigh simulations, the rise in number of 
clumps is, however, much faster when external pressure 
is applied. At later times, the number of clumps flattens 
after ^100 Myr where for the no pressure simulation the 
number of clumps keeps rising with time. 


3.2. Mass Flow Rate 

The star formation rate is sensitive to both the mass 
flux and to the clump mass distribution. We measure 
the gas mass flux through a sphere of radius 16 kpc as 

M = pv -rdS = rrii Vi-Vi ^ , (1) 

iGshell 

where i denotes the index of a cell within a shell of sur¬ 
face S and volume V. In all simulations with enhanced 
pressure, an incoming pressure-driven mass inflow is cre¬ 
ated at the beginning of the simulation and is followed 
by a short mass outflow (seen in the edge-on view of the 
galaxies in Figs. and [^, after which it oscillates around 
zero for the remaining of the simulation (top panels of 
Fig. 0. This is to be expected, as the pressure wave 
coming into the galaxy carries momentum causing the 
galaxy to expel more gas when compared with the non- 
pressurized case where the mass flow rate is close to zero. 

The bottom panels of Fig. shows that the cumula¬ 
tive mass flow remains negative for the gasLow pressur¬ 
ized simulations, indicating a greater mass inflow than 
outflow for these simulations, while the cumulative mass 
flow is close to zero for the non-pressurized simulations. 
Given that the SFR enhancement is larger (factor r\j 10) 
than the net gas change (factor ^ 2), we argue that the 
measured mass inflow does not have a major influence 
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Fig. 1. — Gas density maps (mass-weighted) for the gasLow simulations, for no pressure enhancement (top), and for a pressure enhance¬ 
ment of a factor 3 (bottom), with time evolving from left to right. The galaxies are shown both face-on (upper portion of panels) and 
edge-on (bottom portion of panels). 



Fig. 2. — Same as in Fig. [^for the gasHigh simulations. 
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Fig. 3. — Top: Time evolution of the number of clumps. Bottom: 
SFR as a function of time. The left panel shows the gasLow sim¬ 
ulation set and the right panel shows the gasHigh simulation set. 
Here blue denotes simulations with and green simulations without 
external pressure applied at t = 0. 
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Fig. 4. — Time evolution of the net mass flow rate (top) and cu¬ 
mulative net mass flow relative to the initial mass Mo,init measured 
within a galaxy radius of 12 kpc and a disc height of 1.5 kpc before 
enhancing the pressure (bottom), both at 16 kpc from the galaxy 
center for the gasLow (left), and gasHigh (right) simulations. 


on SFR enhancement measured in the pressurized simu¬ 
lation of gasLow. For the gasHigh simulations, the differ¬ 
ence in the cumulative mass flow is even smaller, show¬ 
ing that the two galaxies have the same amount of gas 
to form stars. 


3.3. Star formation history 

We have seen that an increased pressure enhancement 
leads to an increased number of clumps. Since the star 
formation is proportional to the local density and since 
the gas density threshold of clump detection was set to 
be above that for star formation, one expects a similar 
behavior for the SFR. 

The SFR of the different runs, shown in the bottom 
panel of Fig. for the gasLow (left) and gasHigh (right) 
simulations, resemble the time evolution of the number 
of clumps previously shown in the top panel of Fig. 

In particular, for the gasLow galaxies, while star forma¬ 
tion sets in very late at a modest rate in the simulation 
without external pressure, the SFR rises almost imme¬ 
diately and gradually in the simulation when pressure is 


initially applied. For the gasHigh galaxies, the SFR rises 
gradually starting at early times when no pressure is en¬ 
forced on the galaxy, and later accelerates, similarly to 
the clump behavior. In contrast, when external pressure 
is initially applied to the gasHigh galaxy, the SFR rises 
significantly faster at early times, reaching a plateau af¬ 
ter ^80 Myr, maintaining this rate for the remainder of 
the simulation. 


4. CONCLUSIONS 

We show that a toy model for AGN-induced overpres¬ 
surization leads to enhanced star formation in disk galax¬ 
ies. The effects are dramatic at early times, regardless of 
the initial gas fraction of the galaxy, although the high 
gas fraction galaxy, which is Toomre unstable, experi¬ 
ences a very early rapid burst of star formation, followed 
by more moderate SFR that remains above the SFR in 
the non-pressurized case. However, even in the Toomre 
stable low gas fraction galaxy, external pressure can en¬ 
force fragmentation and star formation, albeit on a longer 
time scale. 

One reason for the increased star formation in pressur¬ 
ized galaxies is early pressure-driven mass inflow from 
the halo and outer disc, on a 10 Myr time-scale (Fig. |^, 
especially on the low gas fraction galaxy. This gas in¬ 
flow feeds cloud/clump growth in the inner parts of the 
galaxy (Fig. ^ and eventually leads to enhanced star 
formation. The infall effects are more dramatic for the 
high gas fraction galaxy, where the disc is more unstable. 
(Fig.I^. We have confirmed this with a detailed analysis 
of the Toomre parameter and clump properties, which 
we defer to a more future study. 

One expects that gas-rich star-forming galaxies should 
have gravitationally-driven star formation rates of the 
order of 100 M 0 yr“^, simply by scaling the gas surface 
density, and maintaining a similar efficiency. Indeed, ob¬ 
servations at z ^ 2 confirm this trend. The more ex¬ 
otic cases of very high (^ 1000 M 0 yr“^) star formation 
rates are more typically found at higher redshift and al¬ 
most invariably have associated luminous AGN. While 
any connection is speculative, and indeed the very di¬ 
rection of possible causality is vigorously debated, our 
simple tests of the effects of AGN-induced pressure, due 
to wind-driven or jet-initiated bow shocks that overpres¬ 
surize the entire inner gaseous disk, suggest that strongly 
enhanced star formation rates are readily achievable. 

Similar SFR enhancements are found from increased 
pressure in the initial stages of ram pressure harassm ent. 


jyx WJ. X lxx xxcxx ctooxxx v^xx u ^ 

although gas loss dominates the long term behavior (|Fu-| 
jit a fc Nagashima 1999 Abadi, Moore & Bower |1999[ 
Bekki||2Q14p . Here, however, the gas is retained, and the 


star formation rate enhancement is far more pronounced. 
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